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Abstract 

In this paper, we develop a variational perturbation (VP) scheme for calculating vacuum expec- 
tation values (VEVs) of local fields in quantum field theories. For a comparatively general scalar 



(N 
> 

in 

en 

field model, the VEV of a comparatively general local field is expanded and truncated at second 



■ order in the VP scheme. The resultant truncated expressions (we call Gaussian smearing formulae) 

in ; 

consist mainly of Gaussian transforms of the local-field function, the model-potential function and 



a* 



their derivatives, and so can be used to skip calculations on path integrals in a concrete theory. 
As an application, the VP expansion series of the VEV of a local exponential field in the sine- 
and sinh-Gordon field theories is truncated and derived up to second order equivalently by directly 



performing the VP scheme, by finishing ordinary integrations in the Gaussian smearing formulae, 
and by borrowing Feynman diagrammatic technique, respectively. Furthermore, the one-order VP 
results of the VEV in the two-dimensional sine- and sinh-Gordon field theories are numerically 
calculated and compared with the exact results conjectured by Lukyanov, Zamolodchikov et al., or 
with the one-order perturbative results obtained by Poghossian. The comparisons provide a strong 
support to the conjectured exact formulae and illustrate non-perturbability of the VP scheme. 
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I. INTRODUCTION 



Green's functions and correlation functions in quantum field theory (QFT), statistical 
mechanics and condensed matter physics are closely related to experimental studies on 
macroscopic matter systems and elementary particles Because of the existence of an 
operator-product-expansion algebra, various multi-point Green's or correlation functions 
with the points approaching each other can be reduced down to one-point functions or 
vacuum expectation values ( VEVs) of local fields [2] . On the other hand, as far as VEVs 
themselves go, they determine the linear responses of statistical-mechanics systems to ex- 
ternal fields, and contain non-perturbative information about QFT, which is not accessible 
through a direct investigation in the conformal perturbation theory. Hence the problem of 
calculating the VEVs of local fields is of fundamental significance. 

There exist some methods for calculating the VEVs of local fields. For an integrable QFT 
which can be considered as a conformal field theory perturbed by some operator, the VEV 
of the perturbing operator can be exactly obtained with the help of the thermodynamic 
Bethe An S at Z approach HQ. This Mad of exact calculation succeeded on ly in very few 
cases J4j. Nevertheless, S. Lukyanov and A. Zamolodchikov made a significant progress 
in 1997 and conjectured the exact VEV of an exponential field in the two-dimensional 
sine-Gordon (sG) and sinh-Gordon (shG) field theories j^J. This progress has led to the 
use of reflection relation [5( and accordingly exact VEVs of local fields in manyperturbed 
conformal field theories without and with boundaries have been proposed through 
solving the reflection relations. Furthermore, these conjectured exact VEVs have again 
been employed to give VEVs of local fields in some other QFTs by making use of quantum 
group restriction existed between the relevant QFTs 0, Q . Besides, recently, because of the 
peculiarity in super-Liouville theories with boundary, conformal bootstrap method |9( and 
Modular transformation method were used to derive one-point functions of bulk and 
boundary operators. For the boundary scaling Lee- Yang model, one-point functions of bulk 
and boundary fields were approximated by using the truncated conformal space approach 
and the form-factor expansion 111 ]. Additionally, in order to check those exact predictions, 
perturbation th eory 12J and some numerical methods based on the truncated conformal 
space approach |l3[ were adopted to calculate the VEVs. 

Obviously, since VEVs of local fields in QFT are non-perturbative objects, general and 
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systematical non-perturbation methods of directly calculating VEVs are needed and such 
tools are also necessary for independently checking those predictions on the exact VEVs. 

n 

We feel that a variational perturbation (VP) theory |l4| can afford such a tool. The VP 
theory can be regarded as some mixture of conventional perturbation theory and variational 
method. As is well known, the perturbation theory is a systematical approximation tool and 
results from it can be improved order by order. Meanwhile, the variational method is feasible 
and effective as well as valid for any coupling (the weak or the strong). These two methods 
have been dominating approximation calculations in theoretical researches for a long time, 
being two standard approximation tools. However, the perturbation theory is valid only for 
very weak coupling at most and the variational method gives no indication of the error in 
its resultant value. In order to collect their merits, avoid their drawbacks, and, of course, to 
develop a systematical non-perturbative tool, a primitive idea of the VP theory for simply 
combining the conventional perturbation theory with the variational method was proposed 
as a tool of solving stationary Schrodinger equation in 1955 [15] (which can even date back 
to even twenty more years earlier ^^). This naive combination has been applied to many 
fields in physics (see references in Ref. Q, Q|). It amounts to an expansion around 
variational approximate result and the variational parameter in the expansion is determined 
with the one-order result by variational method. This primitive VP theory really produces 
non-perturbative results which are valid for any coupling strength and can be improved order 
by order. Nevertheless, the naive combination of the perturbative and variational methods 
does not provide a convergent tool because the variational parameter is independent of the 
approximate order 14]. In 1981 or so, the principle of minimal sensitivity (PMS) [19] was 
proposed by P. M. Steveson (Possibly, earlier in the middle 1970s, a similar principle was 
proposed in Mosc. Univ. Phys. Bull. 31, 10(1976) by V. I. Yukalov) and can be used 
to determine an auxiliary parameter (the aforementioned variational parameter) which is 
artificially introduced into the VP theory. The VP theory with the PMS determines the 
parameter order by order (see next section) and is believed to be a fast convergent theory ^| . 
Now, as a systematical and non-perturbative tool, it has developed with many equivalent 
practical schemes to calculate energies, free energies and effective potentials of systems, and 
has been applied toQFT, condensed matter physics, statistical mechanics, chemical physics, 



and so on [U 



171 Il8j (for a full list, to see references in Refs. [14. 17, [18|]). In this paper, 



we intend to develop a VP scheme to calculate the VEVs of local fields in QFT. 
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In sect. II, we will develop the VP scheme in a general way. For universality and defmite- 
ness, we shall consider a class of systems, scalar field systems or Fermi field systems which 
can be bosonized, with the Lagrangian density [^(j] 

C x = \d^ x d^ x - V{<f> x ) . (1) 

In Eq.(l), the subscript x = (x,t) represents the coordinates in a (D + l)-dimensional 
Minkowski space, dn and are the corresponding covariant derivatives, and <p x = 4>(x) the 



scalar field at x. The potential V(<p x 



in a sense of tempered distributions 



in Eq.(l) is assumed to have a Fourier representation 
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Vfa) = [ -%V{u)e^ . (2) 
J y 2,71 

Speaking roughly, this requires that the integral V(a)e~ Ca2 da with a positive constant C 
is finite. This shouldn't be regarded as a limitation, and, as a matter of fact, quite a number 
of model potentials possess the property, such as the potentials of polynomial models, sG 



and shG models, Bullough-Dodd model, Liouville model, two models discussed in Ref. 22], 
and so on. By the way, a similar general model was studied in Ref. |2j . For comparisons with 
existed work, we will work in a v- dimensional Euclidean space with v = {D + 1). Through 
the time continuation t — > —it with r the Euclidean time, a point x in the Minkowski 
space is transformed into a point r = (x, r) in the Euclidean space. For universality and 
definiteness again, we will also assume that the local field has a Fourier representation, 

Wrb) = f -%6(n)e^° , (3) 

J V 2.71 

at least in a sense of tempered distributions. Here, r*o = (xq, tq) is a given point in the 
Euclidean space. It is evident that such a local field is a comparatively general one. In 
the present paper, based on the VP scheme in Ref. (the scheme stemmed from the 



Okopinska's optimized expansion [23j, and was proposed by Stancu and Stevenson [24]), we 
will develop a VP scheme to calculate the VEV of the local field O((j)f ), Eq.(3), for the 
field theory Eq.(l) in the z/-dimensional Euclidean space. In subsection A of next section, 
we will state the VP scheme. It embraces mainly two key steps: one is the VP expansion 
on the VEVs with an auxiliary parameter introduced, and another is the determination of 
the auxiliary parameter in the expressions truncated from the VP expansion series. It is 
the truncated expressions that give rise to the VP approximate results of VEVs up to the 



truncated order. Then, in subsection B of next section, the truncated expansions of the 
VEV of O(cf)f ) in the field theory, Eq.(l), will concretely be derived up to the second order. 
One will see that the resultant truncated expressions up to the second order are composed 
mainly of Gaussian transforms of the local field 0(4>) and the potential V(4>) as well as their 
derivatives, and we will call them Gaussian smearing formulae. For, at least, any scalar field 
theory which is involved in the model with Eq.(l), one can give VP approximated VEVs 
of O((j)j? ) up to the second order just by finishing ordinary integrations appeared in the 
Gaussian smearing formulae instead of calculating path integrals in the definitions of the 
VEVs. This point is the main reason why we are interested in the general model, Eq.(l) and 
the general local fields, Eq.(3). We should point out that since the QFT with Eq.(l) is not 
a concrete theory, so we will not carry out the other key step to determine the adjustable 
parameter in the Gaussian smearing formulae. Section IV will provide such an example on 
how to determine the auxiliary parameter by the PMS. 

About the VP scheme, the renormalization problem need to be explained here. Since a 
bare field theory is full of divergences and makes no senses, we have to face those divergences 
appeared in the Gaussian smearing formulae and perform a necessary renormalization pro- 
cedure to make the formulae finite before performing the second key step with the PMS in 
the VP scheme. Generally, the renormalization procedure is usually very complicated, and 
is similar to those in perturbative theory [24]. In order to concentrate our attention at de- 
veloping the VP theory, we do not hope to be plagued with the complicated renormalization 
problems of QFT, but we certainly hope to give a finite example for the Gaussian smearing 
formulae. Fortunately, for any two-dimensional scalar field theory with non-derivative in- 
teractions, all ultraviolet divergences can be removed by normal-ordering the Hamiltonian 
2I ■ This fact led to a simple renormalization scheme for two-dimensional field theory, 
the Coleman's normal-ordering renormalization prescription. Furthermore, this convenient 
renormalization prescription has been generalized to path integrals of Minkowski and Eu- 
clidean formalisms in Ref. [3] and Ref. j^J (2002), respectively. Hence, this prescription 
was used in Ref. Q and will be adopted in the present paper. One will see that the Gaussian 
smearing formulae are full of no explicit divergences for the case of v — 2 and so, generally, 
no further renormalization procedures are needed for this case. 

As an application of the scheme, Sections III and IV will calculate the VEV of the local 
exponential field e* ^ in the sG field theory, G a , with a a parameter. The sG model, which 
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appeared early in 1909 27], is closely related to lots of problems in mathematics and physics, 
and has been extensively studied. For the general local field O(<fi? ), Eq.(3), its VEV can 
closely be related to that of an local exponential field, and so the problem of calculating the 
VEVs of the local exponential field is important. In 1997, through direct calculations in the 
sG field theory at the coupling (3 2 — > ( semi-classical limit) and in the free-fermion theory 
(equivalent version of the sG field theory at (3 2 = ~), and through the exact specific free 
energy for the sG field theory, Lukyanov and Zamolodchikov obtained the exact G a for the 
following three special cases: (3 — > (semi-classical limit), (5 = | and a — (3, respectively 4J. 
And then, in the same paper, starting from those exact expressions for the special cases, they 
guessed an exact formula for G a in the two-dimensional sG field theory at any (3 2 < 1 and 
|i?e(a)| < 1/(2/3). Obviously, it is worthwhile checking the conjectured exact formula. In 
order to do so, in the same paper, defining "fully connected" one-point functions, <T2„, from 
the VEVs of even-power fields 2n , they showed that 02 and 04 from the above-mentioned 
exact formula agree with those from perturbation theory for the sG field theory up to (3 A 
and that 02 agrees with the corresponding one-point function from perturbation theory 
up to g the cou plin g in the massive Thirring model, which is the fermion version of the 
sG field theory |25[(1975). Furthermore, it was found that the conjectured exact formula 
can be reobtained from the reflection relations jsj], supporting it indirectly. Slightly later, 
in 2000, checks from perturbation theories in both an angular and a radial quantization 
approaches for the massive Thirring model indicated that the perturbation result up to g 
exaat.y co.nc.des with the corresponding result obtained by expanding the exact fonnn.a 
according to the coupling g 12] (2000). In the same year, a numerical study for the model at 
a finite volume also provides evidence for the conjectured exact formula 13] (2000). In brief, 
up to now, the conjectured exact formula for G a in the two-dimensional sG field theory has 
been completely checked for the case of (3 2 -> \ (g -> is equivalent to (3 2 -> |, See Eq.(67)), 
and received some indirect evidences for its validity. So, besides providing an example of 
applying the VP scheme, investigations on G a will also give the conjectured formula a direct 
check for the cases of (3 2 7^ | (Ref. J^J has provided a partial check for (3 — > 0). 

In section III, we will derive the expressions truncated from the VP expansion series on 
G a at the second order. It will be done in three ways. Subsection A will directly perform 
the VP expansion on G a and then truncate the expansion series at the second order. Since 
the sG and shG field theories are involved in the class of QFTs, Eq.(l), so, in subsection B, 
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the truncated expressions of G a up to second order will be recalculated by finishing those 
ordinary integrations existed in the Gaussian smearing formulae in Section II. The results 
are identical to those obtained in subsection A. Thus, these calculations check and confirm 
the correctness, usefulness and simplicity of the Gaussian smearing formulae in Section 
II. Besides, the VP expansion procedure is formally similar to conventional perturbative 
expansion, and so the Feynman diagrammatic techniques for perturbative theory, which 
have developed very well up to now, can be borrowed into the VP expansion. Subsection C 
will provide such a calculation on G a up to the second order in the VP expansion, and the 
resultant expressions are same as in Subsection A and B. 

The resultant expressions in Sections II and III are not the final VP approximate results 
because the other crucial step of the VP scheme is not yet performed to determine the 
auxiliary parameter by the PMS. As was aforementioned, Sect. IV will do it by considering 
the truncated expressions up to the first order in Sect. Ill and treat only the case of v — 2. 
The sG field theory can be transformed into the shG field theory, and so the one-order VP 
results on the VEVs of the exponential fields in the two-dimensional both sG and shG field 
theories will numerically calculated and be completely compared with the conjectured exact 
results. For the sG model, we will also compare the one-order VP results here with those 
from a perturbation theory in Ref. [12](Poghossian). These comparisons, albeit just the one- 
order results are considered, not only give a strong support for the conjectured exact results 
but also indicate usefulness and non-perturbability of the VP scheme here for calculating 
one-point functions. By the way, the VP approximate result up to the second order on 
the VEV of the exponential fields in the two-dimensional sG field theories has briefly been 



reported in Ref. |28(, and has a less error to the conjectured exact results than the one-order 
result. This gives a sign for the convergency of the VP scheme here. 

Briefly, next section, we will develop the VP scheme by calculating the VEVs of the local 
field, Eq.(3), in the QFT, Eq.(l), and Sect. Ill and IV will provide an application of the 
scheme to the sG and shG field theories. A brief conclusion will be made in Sect. V. 

II. THE VP APPROACH TO VEVS OF LOCAL FIELDS IN QFT 

As was stated in the introduction, we will first state the VP scheme in subsection A, and 
then derive the first three terms in the VP expansion series in subsection B. 
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A. the VP Scheme 



The generating functional is usually the start of discussing a QFT and will be the basis 
for the purpose in the present paper. In this subsection, we will first introduce it, then state 
our VP scheme. 

For any field theory, one can use either the Minkowskian formalism or Euclidean formalism 
to study it 29]. Here, we choose to have our discussion in Euclidean formalism. The 
Euclidean Lagrangian corresponding to the Minkowskian Lagrangian, Eq.(l), is 

C? = ^V r -#V r -0 r - + V{fo) , (4) 

which has the same form with the Hamiltonian density of Eq.(l) in v + 1 spacetime. In 
Euclidean formalism, the corresponding generating functional takes the following form Q|: 

Z[J\ = J £>0exp{- J d v r [~V^V^ + V(fo) - J^f}} , (5) 

where, is the gradient in the ^-dimensional Euclidean space, J? an external source at r, 
and T>(p the functional measure. Generally, for an interacting system, the right side hand 
of Eq.(5) can not exactly be calculated. But for a free field system with V{4>f) = ^fi 2 ^, 
its generating functional is exactly calculable and plays a crucial role in the VP theory and 
relevant calculations. For the convenience of later uses, we write down its result as follows 



a 



= exp{-i J d-fl^jexp^-Jf-'j} . (6) 



Here, 



/ W (Q 2 ) = <| 1 (^) 2 (P 2 +W ' for n^O 

SW? MP 2 + Q 2 )i for n = 



with p momentum and p = \p\ and Jf 1 J = J d u r 1 d v r 1 ' ' J? f^^Jpn with 2^(1990) 



JV zf i ,,vl1—\ 

1 _ / a p 1 _ V 



In Eq.(7), K n (z) is the nth-order modified Bessel function of the second kind. 



S 



Now we address ourselves to the VEVs. For a local field O(cj)f? ), its VEV in the theory 
with Eq.(4) is defined as follows |4| 

In general, O in last equation cannot be exactly calculated and one has to manage to design 
some approximate tool to attack it. 

According to Eqs.(3),(4) and (5), O in Eq.(8) can easily be rewritten as 

= / . (9) 

In the numerator of the integrand in Eq.(9), the subscript means that iQ5(f — rb) takes the 
place of J T in Z[J]. Hence, O can be obtained via calculating the generating functional. 
Luckily, the calculation of effective potential has stimulated the establishment of several VP 
expansion schemes on the generating functional (see Refs. B Q and references therein), 
and they can possibly be developed to calculate O. Next, we will g eneralize the VP scheme 
in the Ref. 18] (a very slightly different version of that in Ref. |24j ) to calculate O. 

As was stated in the introduction, the Coleman's normal-ordering renormalization pre- 
scription will be adopted. That is, the Euclidean Lagrangian in the exponential in Eq.(5) 
will be replaced by the following normal-ordered form with respect to an arbitrary normal- 
ordering mass M Q (2002), 

-\k0)(M 2 ) + ~M% } (M 2 ) + ~Vr0rV r -0 r -- J r ^+ [ -^v(uy^ e ^k^) 
Z Z Z J \jZix 

In the above expression, the last term is the normal-ordered form of V(</y), and can be got 

by using Eq.(2) and the Baker-Hausdorff formula e A+B = e A e B e~^ A ' B ^ with the commutator 

[A, B] some c-number. So, Z[J] is turned into the following form : 

Z[J] = exp{ J d»r[ll (0) (M 2 )-±M 2 I {1) (M 2 )]} J V<P 

xexp{- / d v r [IvrfrVrt?- [ ^V(oo)e iuJ ^e^ 2l ^ M2) }} . (10) 

J 2 J w2n 



Eq.(10) is nothing but Eq.(ll) of Ref. 20] (2002) in notations here. Thus, for the case of two 
dimensions, the fields and the model parameters, such as mass and couplings, are now finite 



(for simplicity, we use the same symbols as those in the bare 



Now we further modify Z\J\ in Eq.(10) by following Ref. Q] or 



-agrangian, Eq.(4)). 
2o| (2002) (only without 



shifting the field, for simplicity). First, a parameter \x is introduced by adding a vanishing 



term J d u f^(j)ff(/j, 2 — fx 2 )(f>f into the exponent of the functional integral in Eq.(lO). This way of 
introducing the artificial /i was used several decades ago jscj . Then, rearrange the exponent 
into a free-field part (with a mass being the introduced parameter fx) plus a new interacting 
part. Finally insert a formal expansion factor e in front of the interacting part. Consequently, 
Z[J] is turned into the following form 

Z[J;c] = exp{|rfV[^ (0) [M 2 ]-^ 2 / (1) [M 2 ]]} 

x J P0exp{- J d u f^M-V? + » 2 )<f>?+ eni{<f>e, - J r fc]} 
= exp{ J d"r [^/ ( o)[M 2 ] - ~M 2 /(i)[M 2 ]]} 

x J P0exp{-e J dVWj(— ,/i)}exp{- ^ d"r[^(-V£ + // 2 )<fc - J r 4r]} 
= exp{- | drr[^ ) (j^-^(M 2 ) + ^M% ) (M 2 )^ 

xexp{-e / ^fW J (A ;/ ,)} exp {Ij/- 1 J} (11) 
with the new interacting part 

= -^ 2 #+ / -^VC^e^ei^w^ . (12) 

^ J V 27T 

In writing down Eq.(ll), we have employed the result Eq.(6). Note that extrapolating Z[J; e] 
to e = 1, one recovers Z[J] in Eq.(lO). 

Expanding exp{— e J d v r 7~ti{jjz, A*)} i n Eq.(ll) into a series in e, one has 

Z[J;e] = exp{- J d»r ~ ^(M 2 ) + \m 2 I {1) {M 2 )}} 

x E en Hl! / Yl^M^rieM-Jf-'j} . (13) 

^ — n ' ^ / — 1 r k 



Thus, we obtain an expansion of the generating functional without requiring the model cou- 
pling and/or the formal expansion factor e very weak. It seems that this expansion form 
has no senses and is useless. Nevertheless, when the PMS is introduced, it can really fur- 



ms 



i quite a potential non-perturbative approximate tool for us to attack physical problems 
14, Q, Q]. Next, one will see that it provides a basis for non-perturbatively producing 
approximate results of O. 

Since we are using normal-ordered form, the local field C(0r o ) should be replaced by its 
normal-ordered form from now on, and, from Eq.(3) as well as the Baker-Hausdorff formula, 

10 



it can be written as 



V Z7T ^ 



at least in a sense of tempered distributions. Then O in Eq.(9) can take the following form 



O 



^0(n)ex V {-n 2 I {1) (M 2 )} 



/ flLi W/djr. A*) exp{i J/- 1 J}] J-= 



x- 



in5(f— fo) 



J=0 



(14) 



e=l 



The right hand side of Eq.(14) is only a quotient between two series in terms of e and can 
be rearranged into a series consisting of powers of e 



O = [C (0) + e l O {1) + e 2 (2) + ■ ■ ■ + e n O {n) + ■ • ■ ] e=1 = I e n O (n 

n=0 

This can be realized according to the formula 0.313 on page 14 in Ref. 

^ oo 1 n 

^ C-nZ ) C n -|- ^ C n —kdk b n - 



e=l 



(15) 







«o 



«o 



=0 ttnZ ' -« n =0 -» fc=l 

This series, Eq.(15), is the VP expansion series with the auxiliary parameter fi. For v = 2, 
every terms in the series Eq.(15) must be full of no explicit divergences (there are possibly 
some implied divergences for some theories, and such an example will be met in next section), 
and truncating it at some order of e will lead to an approximate O. In the case of v > 2 
where divergences are met in Eq.(15), the model parameters (the mass, coupling constants) 
should be written as series in e as did in Ref. , respectively, and then one can successively 
do as follow: substitute the series of the model parameters into the truncated expression, 
rearrange the truncated expression in terms of powers in e, keep terms in it only up to the 
truncated order. After these, one can take e = 1 in the truncated expression and manage 
to find a renormalization scheme for rendering the truncated expression explicitly finite 
(finding the renormalization scheme is similar to what is done in a perturbative theory). The 
finite truncated expression at any order will be dependent upon fi. How can we determine 
the arbitrary parameter [i in a truncated result? As was stated in the introduction, we 
can determine /i with the PMS jl^. Evidently, the exact result, Eq.(15) with e = 1, is 
independent of the auxiliary parameter /x and is a constant in the space of /i. Hence, it 
would be reasonable that the truncated result should vary most slowly with the parameter 
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/i so that it can likely provide a most reliable approximate result to the exact one. This is the 
main spirit of the so called PMS A simple realization of the PMS is to require the first 
derivative of the truncated expression with respect to /i to be zero. Sometimes, it cannot 
give rise to a meaningful solution, and in this case one can render the second derivative zero 
to determine fi. Thus the above procedure, truncating the series in Eq.(15) at some order 
of e and using the PMS to determine /z, provides an approximate method of calculating O 
which can systematically control its approximate accuracy. 

In a general, we can say nothing on the convergent property of Eq.(15), let alone on 
the convergency of the sequence consisting of the truncated results at various orders, be- 
cause it is difficult to prove their convergency in a general way. Nevertheless, from the 
aforementioned spirit of the PMS, the convergency of the sequence of the truncated results 
determined with the PMS to the exact result should be conceivable and understandable. 
Existed investigations and applications have indicated and illustrated the point. In the 
early 1980s, Stevenson proposed the PMS, studied and indicated the efficiency of the PMS 
for keeping the VP theory convergent. In 1990s, it was rigorously proved that it is owing 
to the PMS that the VP theory leads to quickly convergent results on ground state energy 



14j(2004), and recently, the rigorous 



for a quantum-mechanical anharmonic oscillator 
proofs of the convergency on the VP theory have been given to the critical O(N) scalar field 
theory j^. From these existed investigations, the VP scheme of calculating VEVs in the 
present paper is presumably convergent, because, frankly, it is just a generalization of the 
existed VP schemes of calculating energy, free energy, effective potentials, and so on. As 
was aforementioned, a brief report in Ref. [28] suggest a sign on the convergency of the VPT 
scheme here. In the present paper, we will not discuss the convergent problem. 

The VP scheme of calculating VEVs stated in the above is a non-perturbative method. 
The expansion in Eq.(15) in terms of e is intrinsically different in nature from a conven- 
tional perturbation expansion in terms of the interaction coupling, albeit they are formally 
similar to each other. The key point of the reliable approximate results from conventional 
perturbation expansion consists in the requirement of weak couplings, whereas the crucial 
point of the reliable and fast convergent approximate results from the VP expansion consists 
in the requirement of most insensitivity of the truncated results to the auxiliary parame- 
ter. It is due to the adjustability of the auxiliary parameter that the expansion in Eq.(15) 



has no limitations both to the formal expansion factor e (in fact, e = 1 is taken eventu- 
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ally) and to the coupling. In the VP expansion, there exists the adjustable term, — |/i 2 2 , 
besides the original interaction terms (see Eq.(12)), and the propagator Eq.(7), which is 
embraced in the unexpanded part of Z[J; e], also contains the adjustable parameter, \i. So 
it is conceivable that for a different value of the coupling (weak or strong), one can adjust 
the value of n through the PMS to get to a reliable approximate result of O. The existed 
investigations have indicated the non-perturbability of the VP theory. In Section IV, we 
will provide a comparison of the VP theory here with perturbation theory to illustrate the 
non-perturbative property of the VP scheme here. 



B. Truncate the VP Expansion Series at the Second Order 



This subsection derives the first three terms in the series Eq.(15). At the zeroth order, 
Eqs.(14) and (15) lead to 



0(0) 



exp{-J/-V} 



Jp=iQ5(r— ro) 



2tt 

J-oo V27T 



2tt 



00 da e _a? +inaVU 



2tt 



(16) 



with U = J( 1 )(/i 2 ) — I(i)(M. 2 ). Here, we have used Eq.(3) and the Gaussian integral formula 
-^ e ~^" +av ^ 2a = e a , which will repeatedly be used for obtaining next Eqs.(20) and (25). 
At the first order, the coefficient O^- 1 ' is as follows 



^d(fi)exp{ifi 2 / (1) (Al 2 )} 



J rfV 1 H / (^ 7 ^)exp{ij/- 1 J} 



+ 0(0) 



Jft=iilS(r— ro) 

(17) 



.7^=0 



To calculate {l \ one can first have 



^-exp^JZ-V} = [/"}. + (J/- 1 ) 2 ] exp{ijrV} 



ri 



and 



5 1 1 

exp{^^j-}exp{-J/ _1 J} = exp{-J/ _1 J} 



Jp+iu)5{f—ri ) 



(18) 



(19) 
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with (J/ -1 )™ = (/ d v rJ?fzlJ m . Eq.(19) is easily obtained by returning its left hand side 
to its original integral expression. (So are Eqs.(23) and (24)). Then substituting them into 
Eq.(17), one can obtain 



2tt 

/dVL du f 
V2n V2n J 



d v nO(n)V(u)e— 
duo 



-J 1 



da 



d u nV(u)e-^ 



-oo v27T 



e 2 



^(tVw) 



x / ~7=e 2 

oo v27T 



J rnr 



where C (n) (; 



d n O(z) 



f -^=d{Q)(iQ) n e inz and W jfe = - 



Finally, Eqs.(14) and (15) give the second-order coefficient in Eq.(15) as 



0(2) = 



1 



1 



o(n)ex P {-n 2 i {1) (M 2 )}- / d'r^mA^r^) 



2! 

xWj^^exp^J/- 1 ,/} 



Jp=iQS(r— ro) 



-±<9<°> 

2! 



J ? =0 



+<5< 1 >[ / «^(-^,/x)exp{ij/- 1 J} 
Similarly to what was done at the first order, one can first have 



(20) 



(21) 



1^2 

2 



+/^ 2 (^r 1 )i + (^r 1 )i(^r 1 )2}ex P {-j/- i j} , (22) 



exp^^-}exp{ij/-V} = {[/^ + (J/^)?]exp{ijrV}} 



Jp— » Jp+iLoS (f—r<2. ) 



(23) 



and 



(5 (5 1 1 

exp{^i— -}exp{^2— -}exp{-J/ _1 J} = exp{-Jf^J} 



8 J? 



' 5Jf? 2 ' L 2 



icji <5 (r — r 1 ) + iu>2 <5 (r — fb ) 



, (24) 
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and then substituting Eqs.(22), (23) and (24) into Eq.(21) leads to 

+ 2^ /r - i/ -Je-^- 2 )^-^ 2 
1 f „,_ /" do;2 ^./^x-, , - / 

+ 2j d r ^J W*W*VS°^ VMVM 



' r\ri> 



-^/^[2(i:y 2 + (^) 2 ] 

+£<5<°> / Wf 2 | ^^M[/^ i + M 2 (/-) 2 ]e -^ 

-law [d^f 2 /^^(^(^e-^^^ 
2 J J v^v^ 

r 1 



^ 2 



r-1 

x 0(a v ^ + 7^)-0(av / W)] 

V / drW^f^ £ V«( 7 v^)0«(a>/V^ + 7^ 

- l -!d^r 2 r d -^pBe-^-^V{X^)V{ 1 ^ 2 + A)\0{aVU) 

2 7 J-oo V87T 3 L 



Wo2 ^ + 7 wv^ W 



X 
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f-i , 
O(a v / VHu" + 7^)-0(«v / W) .(25) 



Notice that in Eqs.(20) and (25), those expressions in all the parentheses which follow the 
symbols V, V^ n \ O and are arguments of the functions 0(z) and 0^ n \z), respectively. 
From Eqs.(16), (20) and (25), we can write down truncated VP expressions of O at the first 
order, O 1 = + eO^ , and at the second order, O n = + e(5« + e 2 ( - 2 \ 

It seems that Eqs.(20) and (25) (so O 1 and O u ) have space-volume divergences because 
both O(aVU) and 0^ 2 \aVU) are independent of space coordinates fi and r 2 . However, 
one can have a simple analysis and see that the divergences really cancel. Formally, one can 
write down the following expansion form 

f— 1 oo 

CW^ + T^) = 0(aVU) +Y.<& > 

* n=l 

where the coefficients functions of a, 7 and U but independent of the coordi- 

nate fi, and evidently, the similar expansion forms one can have for (2 \a^W^ 2 + 
7 ^f),0(av/yV^ + 7 ^),0( 1 )(aVH^ + as well as V^jW^ + X%-). Besides, 

one can also has 

/ (Uf -1 — f _1 f _1 ) 2 lAf~ l — f _1 f _1 f _1 

00 

i,m,n=0 

where l,m,n in the summation are required not to take zero simultaneously, and the co- 
efficients a; )mj „s are independent of the coordinates fi and f* 2 . Thus, there exist really no 
space- volume divergences in Eqs.(20) and (25), respectively. Consequently, both O 1 and 
O n ) have no space-volume divergences. 

Furthermore, substituting the above formal expansion forms into Eqs.(20) and (25), one 
will find that in Eqs.(20) and (25), all the integrals over fi and f 2 are involved only in the 
following types 

n-l 

dy,dy 2 ■ ■ ■ dy n ^[(2n)^( P l + ^){ P \ + ^ • • • ^ + //)(( J>) 2 + ^ » 
i (l ' m ' n) = J ' ^n^W^Yif^rif^r = jdy 1 dy 2 ---dy l dy 1 dy 2 ---dy m _ 1 

l m—l 

xdy\dyi 2 ■ ■ ■ dy> n ^) [l+m+n - 2> + ^ + ^ 

i=i 3=1 
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n—1 I m—l / n—1 

fc=l 2=1 j=l 1=1 fc = l 

and products of 7^s. In the case of v = 2, W is finite, /C' m ' n )s as well as with any 
i, m and n are finite, and accordingly Eqs.(16), (20) and (25) are really full of no explicit 
divergences. Since both O 1 and O 11 are finite when v = 2, we can take e = 1 and perform 
the other key step of the VP scheme to termine the parameter fi with the PMS. For the 
case of v = 2, according to the PMS, rendering vanishing, one can solve it for /i = /i 7 
which is the value of /i at the first order, and rendering vanishing, one can solve it for 
/i = fi 11 which is the value of \i up to the second order (in case the first derivative condition 
couldn't produce meaningful root, the second derivative condition T^jp — will be used). 
Substituting fi 1 into O 1 gives the approximate result of O up to the first order in the VP 
scheme, and substituting \x n into O 11 will give the approximate result of O up to the second 
order in the VP scheme. Ref. j3] gave an example for doing so. 

Presumably, although there are no explicit divergences in the case of v — 2, for non- 
polynomial potential V(4> x ) and/or non-polynomial-type local fields O(4>^ ), there possibly 
exist some terms which is not convergent for some range of model parameters in the approx- 
imate results of O. Furthermore, for v = 3, although U is still finite, but 1^ with any n > 2 
and /( i ' m ' n )s with / + m + n > 5 are no longer finite, and for v = 4, both U and J*™** with 
any n > 1 and /C' m ' n )s with / + m + n > 3 are divergent. For all the divergent cases, one 
have to appeal to a further renormalization procedure before determining [i with the PMS, 
as we have stated in last subsection. In the present paper, we do not discuss it concretely. 

In the same way, one can consider higher-order cases. Here we do not continue it. In 
this section, we have developed a VP scheme of calculating VEVs in QFT, and obtained 
the truncated VP expressions for O at the second order, Eqs.(16), (20) and (25). The 
right hand sides of them are mostly Gaussian transforms of the functions O(z), V(z) and 
their derivatives |34j, and Eqs.(16), (20) and (25) are the aforementioned Gaussian smearing 
formulae. According to these formulae, one can easily obtain truncated VP expressions for 
VEVs of a local field in a field theory by finishing only ordinary integrations. 

III. SG FIELD THEORY: TRUNCATIONS FROM THE VP SERIES 

This section will derive the first three terms in Eq.(15) for G a in three ways. 
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A. Direct Performing the VP Expansion 



We consider the //-dimensional Euclidean sG field theory with the following Lagrangian 
density 

£sG — -Vr0rV^0^— 2Q cos(v / 87t/3</v) . (26) 

The Lagrangian density Eq.(26) is nothing but Eq.(5) in Ref. |4| if one makes the transform 
— > (hereafter, we will use e ^ ^/8W , ( ^ o) instead of e ia ^ r °^ as the form of the local exponen- 
tial field and consequently the parameter a and the coupling (3 in this paper are identical to 
those in Ref. Q], respectively). If taking y/87r@ — > [3 and 2Q = m 2 / (3 2 and adding the term 
m 2 / (3 2 in the Lagrangian density, one will get the Euclidean version of the sG Lagrangian 
density which discussed in Ref. |l8j . Besides, takin g th e substitution (3 = and — > —A, 
Eq.(26) describes the Euclidean shG field theory p,|35j|, and so the resultant expressions in 
this section can easily be used to give the corresponding expressions for the shG field theory. 

In Eq.(26), (3 is the coupling parameter with the dimension [length] and Vt is 

another parameter with the dimension [length] _jD_1 in natural unit system. It is always 
viable to have (3 > without loss of generality. Obviously, the classical potential V{4>r) = 

2-wn 



— 2VI cos(y/8ir (3(f) is invariant under the transform 0^0+ with any integer n, and 



so the classical vacua are infinitely degenerate. So do the quantum vacua, as was shown, 
for example, by the beyond- Gaussian effective potential for two-dimensional case Q|. Here, 
we choose to consider the vacuum with the expectation value of the sG field operator </y 
vanishing. 

According to the definition Eq.(8), the VEV of the local exponential field e ly ^ a ^ ro ' > in 
the sG field theory is defined as follows 

G =< e lV ^ am >= / V< ^ exp{iy / 8¥a0(O)} exp{- / d u rC sG } , . 

J Vcj)exp{- J d»rC sG } ' 1 ' 

For simplicity, the exponential field in Eq.(27) is taken at r*o = 0. It is evident that the 
numerator and denominator in the right hand side of Eq.(27) can be easily got from the sG 
generating functional 

Z sG [J] = J V<j)exp{- J d v f[C sG - Jr4?]} , (28) 
by taking J? = i\f&na5{f) and J = 0, respectively. 
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Doing as was done for getting to Eqs.(14) and (15) in last section, one can obtain a 
similar expansion for G a . In the present case, all the things are same as in subsection A 
of Sect.II, except for, now, M M [O{4>r )] = exp{47ra 2 / (1) (M 2 )}e^ a ^ ) and M M \V (<!>?)] = 
— 2f2exp{47r/3 2 /(i)(A4 2 )}cosv / 87r/30^. Consequently, one can have 



K~ o ^ I YIU d»Wf{^ n) exp{| J/- 1 J}] 



Q = e 4na 2 I (1) (M 2 ) 



/87ra5(r) 



[EZo ^ I ECU d»Wf(d-, n) exp{| J/-i J}] J=Q 



e=l 



with 



ftf (<fo /x) = -^ 2 # - 2fi cos(v / 8vr/30 r -) exp{47r/5 2 J (1) (A4 2 )} 



and the series consisting of powers of e 

oo 

G a = [Gf + e 1 ^ 1 ) + e 2 Gi 2 ) + • • • + e »Gl») + • • • W = £ "GW] 



(29) 
(30) 

(31) 



n=0 



Next, we will derive the first three terms in Eq.(31). 
At the zeroth order, Eqs.(29) and (31) lead to 

Gf = exp{4vra 2 / (1) (^ 2 )}[exp{ij/- 1 J} 



- Jp=i-^/8naS(f) 



= e 



-47ra 2 W 



(32) 



At the first order, the coefficient G^ is as follows 

G W = - exp{4 na 2 I {1) (M 2 )}\ [ d^f 1 Hf(^,fi)exp{l-Jf- 1 J} 



+G<°> 



To calculate one can first have 
^f(^,/i)exp{ijrV} 



5J ?1 

sG ^,/i)exp{ijrV} 



J?=i\/8Tra5 (f) 



(33) 



1 



1 



-2^^exp{-J/-V} - fiexp{47r/3 2 / (1) (A4 2 )} exp{- J/" 1 ^! 



ri 



I J^+i \/8nf36 (r — r*i ) 



+ exp{-J/- 1 J}|. 



• Jp— i \/87r/3<5 (r — r*i ) 



(34) 



In writing down Eq.(34), we first returned the cosine part of its left hand side into the 
original functional form, and then used exponential form of cosine function and the result 
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Eq.(6). Substituting Eq.(34) together with Eq.(18) into (33), one obtains 



^^)T 



/ «(/o>I) 2n 



(35) 



Finally, we consider the second order. From Eqs.(29) and (31) , one can write down the 
second-order coefficient in Eq.(31) as 



/7( 2 ) - 1. 47ra 2 / (1) (.M 2 ) 

a ~ 2! 



5.1 



n 



'T1 



d»nd»mf(^,fi)Hf(j^,fi)exp{±Jj-,n 



Jf=i\/STraS (f) 



2! a 



fr^mf (t^, ^Wf (t^, /i) exp{i J/- 1 J} 
J^ 2 z 



J J P =0 



+G«[y «?Yf( 7 ^,/i)exp{ij/- 1 J} 



J F =0 



Doing as was done at the first order, one can first have 

Hf(-£-> m?{^r, U) expi^J/- 1 J} 
o J;?, o z 



ri wu r2 

5 2 5 2 1 



4- 5J3 



ynexp{47r/3 2 / (1) (7W 2 )} 



5Jl 

n 



exp{-J/-V}|. 



1 



+fi 2 exp{87r/? 2 / (1) (A^ 2 )} 



+ exp{-J/- 1 J} 

e^djf-^jyL 



+2exp{ijr 1 J}| 



Jp^Jf!—iV8nPS(r—r2) 



J ? +iV87r/35(r-ri)+iV87r/35(r-r2) 



J?^ J?+i (r - r 1 ) - i v / 8^/35(r -r 2 ) 



+ exp{-J/- 1 J}L 



• J^— iV8n{35(r— ri)— i\/8nf36(r— r*2) 



Substituting Eq.(22) into Eq.(37) leads to 



^^-,/.)Hf(^-, /U )ex P {ij/-V} 



(36) 



(37) 



+^{jr 1 wr 1 )2 + (j/- 1 )?*.//- 1 )!} ex P {i j/- 1 j} 

+ / i 2 fiexp{4vr/5 2 / (1) (^ 2 )}{2[/^ 1 + (J/ -1 ) 2 + 8vr(^) 2 (/ r ^) 2 ] cos[V&F/?( J/- 1 ) 2 ] 

-4v^F)9/^( J/" 1 )! sin[v / 8^/?(J/- 1 ) 2 ]} exp{i J/- 1 J + 47r(^) 2 (/ r ^ 2 ) 2 } 
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+n 2 exp{8rrf3 2 I {1) (M 2 )} 
1 

2" 

+2 cosIV^G//- 1 )! + (J/- 1 ) a )] exp{i J/" 1 J + 87r(0) a [/^ + }} • (38) 



x{2exp{-Jr 1 J + zv / 8T/3[(Jr 1 ) 1 - (J/" 1 ) 2 ] +87r(z/3) 2 [/ r ^ - fz\,}} 

.1 
2' 

For simplicity, in writing down last equation, we have used the the property that Eq.(36) is 
invariant under interchanging r\ and r 2 - (By the way, for getting Eqs.(13),(14) and (T6) in 
Ref. [ill, we also used the similar property, and in the right hand side of Eq.(16) in Ref. jlq . 
the first term should have an additional factor " | e J f lj//2 " and the second term should 
have an additional factor " ~ ".) Thus, substituting Eqs.(32),(35) and (38) into Eq.(36), 
one can eventually obtain 

-&<K f i 2 na 2 e-^ a2 +^ u J (fnij-^f J rfV 2 [cosh(87ra/3/ ^)-l] 
-8n^ 2 Qp 2 e~^ a2+ ^ u J <Fr x d»r 2 {f^ ) 2 [cosh(87ra/3/-i) - 1] 
+lQn^ 2 na(3e- 4 ^ a2+ ^ u J d»r x d v r 2 f^J-\ sirm(87ra/?/-i) 
+2fi 2 e- 47r(a2+2/32) "| y ^fi[cosh(87ra^/--.;) - 1]}' 

+n 2 e-^ a2+2 ^ u J d v r x d 2 r 2 [e~^ r ^ - 1] [cosh(87ra/3(/-J + f w \)) - 1] 
+n 2 e -Ma 2 +2^)u j d ^ lrf v 2 [e 8 ^ 2/ ^2 - i][ e - 8 ^(/o>i-/o-4) _ i] . (39) 

From Eqs.(32), (35) and (39), one can easily write down = G^ + eGa and G 1 / = 
G^ + eG^+e 2 G^\ which are the truncated expressions from the VP expansion series, 
Eq.(31), at the first order and the second order, respectively. It is evident that Eqs.(32), 
(35) and (39) with v = 2 have not any explicit divergences, providing an example for our 
analysis in last section. Using the result U = — ^Mn(-^p-), Eqs.(32), (35) and (39) with 
v = 2 can lead to Eq.(12) in Ref. Q- 

B. Using Gaussian Smearing Formulae 

This subsection will substitute the concrete expressions of the sG potential and the lo- 
cal exponential field into Eqs.(16), (20) and (25), respectively, and finish those Gaussian 
transforms to recover Eqs.(32), (35) and (39). 
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First, from Eq.(16), G^ can be written as 



G (0) = / _^ e -^ e iV8taaVU (4Q) 



Employing the Gaussian integral formula between Eqs.(16) and (17), one can easily check 
that last equation really gives the result Eq.(32). 

Second, according to Eq.(20), takes the following form 

J J-oo V27T 



da _<£_ n 



x / ^=e"^ 

-oo v27T 



(41) 



Finishing integrations over a and 7 in last equation with the help of the Gaussian integral 
formula, one can have 

G« = -47r/iV y d v r x {j-}^e-^ u + 2VlJ d u f ie -^ a2+ ^ u [cosh(87raf3f }) - 1} . (42) 

Expanding cosh(87ra/3/^) in Eq.(42) as power series leads to nothing but Eq.(35). 

Eq.(39) can also be obtained from the Gaussian smearing formulae. Substituting V(<j> x ) 
and O((f)f 0= o) for the sG field theory into Eq.(25), one has 

da a 2 



+8 7 rV^Y««(/^ 1 ) 2 (/^ 2 ) 2 / 



00 V27T 

00 j 2 



■DC' 



2tt 



fZ 1 *, 



x |gi«v / 8^Atea ^iv / 8^(a^+/3v / W)7 + giv / 8^(«^-/3v / W)7j 
_giav'87rWa |~gi/V87rW7 _|_ g-j/V87rW7j | 

J ' ^ J-ooV^V^ 

/"I /"I 

x jgiav^HAtea ^v^(a^3.+M7) 7 + gi> /fcF(a-^-/MJ)7] 

_ e iaVSTvUa ^ipVEnU-y _|_ g-j/V87rW7j j 



2 2 

27T V27T " 
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X( ,iaV*nW<Y2n y v — w , .... _ ( "-^77 



da d'j dX 



_|_g «,/ w 12 



ia\/8irUa 



-fr / d v r\d v r 2 \ —^=—^=—^=e 2 2 2 



2 



2vr V 7 ^ V27T 

/-I /-l 

_ e «av / 87rWa ^^v^ttM _|_ e -i7/V87rWj ^iA^v 7 !^ _|_ g - iA/3 v / 8^i7j j 



Carrying out all integrations over a, 7 and A in last equation with the aid of the Gaussian 
integral formula, one can get to 



7i ) 2 

r r 2 / 



-SnaVShT 4 ***?* J d v nd^r 2 U7Xf[^M^^X) ~ ^ 
-87rp 2 ^Qe-^ a2+ ^ u j d v r,d v r 2 {fzl) 2 [cosh(87ra/3/ry - 1] 
+ 167ra/?/i 2 fie- 4 ^ a2+ ' 32 ) w | d v r x d v r^JZ^ sinh(87ra/3/ry 
^V 4 ^ 2 ^ | cosh^ 2 /^) 

_ 4fi 2 e -4.(a2 + 2 /3 2) W y ^f 1 ^f 2 [ cosh ( 87ra/ 5/-l i ) _ 1] 

: <k ) j . (44) 
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Noting that the coordinates fi and r*2 have equivalent positions, one can easily find that 
Eq.(44) is identical to Eq.(39). Thus, instead of directly calculating path integrals, we have 
reobtained Eqs.(32), (35) and (39) only by finishing ordinary integrations in the Gaussian 
smearing formulae. 



C. Borrowing Feynman Diagrammatic Technique 

In this subsection, we show that the results Eqs.(32), (35) and (39) can be reobtained by 
borrowing Feynman diagrammatic technique with the propagator Eq.(7). It is well known 
that one-point functions in the perturbation theory are the sum of all connected Feynman 



diagrams 36[. This point is also true for Eqs.(15) and (31) because the VP expansion is 
formally similar to the perturbative expansion, and accordingly in Eq.(31) is only the 
sum of all nth-order connected Feynman diagrams constructed with an external vertex aris- 
ing from exp{iv / 87ra0(O)} and n internal vertices arising from 7-^f G (</y, «). The Feynman 
diagrammatic technique for the sG perturbation field theory has developed early in 1970s. 
For example, to perturbatively check the equivalence between the Coulomb gas and the sG 
field theory, Feynman diagrammatic technique was established in 1978 [37], and Feynman 
diagrams in the sG perturbation field theory up to the third order were also analyzed to con- 
sider the renormalization problem in the sG perturbation theory [38] . In our VP expansion, 
the Feynman diagrammatic technique is similar to the perturbation theory, but, differently, 
the propagator is Eq.(7), and vertices are classified as external vertices which come from the 
exponential fields and internal vertices which come from Hfifa/J,) in Eq.(30). 

Expanding the exponential field e lv ^ a ^ in powers of 0(0), one will find that there 
exist infinitely many external vertices (we will call them e n - vertices), having any number 
of legs, and the coefficient e n = e 47 ™ 2 ^!)^ 2 ) (gvgjp) adhere to the external n-leg vertex, 
the e n -vertex. Furthermore, Expanding the cosine field cos(v / 8~7r/30r) of Eq.(30) in powers 
of 4>(f), one can have infinitely many internal vertices (called Cm- vertices here), having 
any even number of legs, and an internal 2n-leg vertex, C^n-vertex, possesses the coefficient 



C 2n = -2ne^ 2l w (M ^ {iV ^f n . Besides, the first term of Eq.(30) contributes an additional 
internal 2- leg vertex (called u- vertex) which the coefficient U2 = — |u 2 adheres to. With all 
the vertices and the propagator, one can draw out all the Feynman diagrams order by order 
for G a and obtain a diagrammatic expansion of G a which is identical to Eq.(31). 
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The zeroth order diagrams for G a are all possible connected graphs self-contracted by the 
e n - vertices. The external vertices with (2n+l) legs cannot be completely self-contracted and 
have no contributions to G a , which correspond to the vanishing contributions from the odd 
powers of the fields. Thus G a °^ is the sum of all possible connected graphs self-contracted 
by the eon-yertices. We draw them in Fig.l (Feynman diagrams are drawn with JaxoDraw 
package [39|] ) . In Fig.l, an n-loop graph has a factor ^\^ n which is the number of ways to 



a 

eo 



+ 



+ 




+ 



n loops 



+ ...= <g> 



64 



e 2 



n 



FIG. 1: The Feynman diagrammatic form of G a 



(0) 



self-contracted legs of the e 2n -vertex. From this graphic form, one can have 
G a -2^ e <2n)\ Uo ° j X , 

n=0 

Obviously, Eq.(45) is identical to Eq.(32). 



n\(2\) n 



(45) 



Now we consider the first order diagrams for G a . They are all possible connected graphs 
each of which consists of any e n -vertex and any internal vertex. The diagrams formed by an 
external (2n + l)-leg vertex and an internal vertex are not completely contracted because 
any internal vertex has even number of legs, and accordingly Ga is the sum of all connected 
graphs each of which is constructed by an e 2n -vertex and an internal vertex. These graphs 
have two types: each of one type consists of an e 2n -vertex and the /z-vertex (e 2n -Ai graph) 
and each of another type consists of an e 2n - vertex and a C 2n - vertex (e 2n -C2n graph). We 
first consider the sum of the e 2n -/i graphs. It is shown in Fig.2. In Fig. 2, an e 2n -/i graph has 



e 2 ^2 e 4 112 [Jcq 112 





(n — l)loops 



+ ••• + ' 




FIG. 2: The sum of all connected first-order graphs consisting of e 2n -vertices and the /i-vertex. 
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a factor ^ n _ 1 j-, 2 ^ ! (n _ 1) . The sum in Fig. 2, G^\, has the following result 

(1) _ r-l\n-l„4Tra 2 I a) (M 2 ) 0^/8^0^ f , v -*(f-\\Z( 1 ,,2\ ( 2n )' 



(n- 1)!(2!)™" 1 

= ^(iV^af J d^f^fG^ . (46) 

Note that at the first order of e, there is an additional total factor ^— p- the expansion factor, 
and we have added it to Eq.(46). We do so in the later Eq.(47), too. 

Comparing the result Eq.(46) with Fig. 2, one can see that the sum of all possible con- 
nected graphs consisting of the e2 n - vertices and the u-vertex amounts to the sum of all graphs 



consisting of an external 2-leg vertex with a coefficient ( iv ^ a ) q@) (fa e vertex represented 
by the cross-centered circle with two legs in the right hand side of Fig. 2) and the «- vertex 
(in fact, only two of such graphs). Noting that the //-vertex has two legs, one would admit 
that the sum of all possible connected graphs consisting of external vertices and an n-leg 
internal vertex amounts to the sum of all graphs consisting of an external n-leg vertex with 
a coefficient (jV ^ a) " Gi 0) (we will call it an external n-leg exponential vertex) and the n-leg 
internal vertex (one can prove it directly in the way of getting to Eq.(46)). Similarly to 
those in Fig.l and Fig.2, we will graphically represent an external n-leg exponential vertex 
by a cross-centered circle with n legs. Furthermore, one can also check that the sum of all 
connected graphs which consist of C^-vertices and some external or internal 2n-leg vertex 
amounts to summing up all connected graphs which are made by an internal 2n-leg vertex 
with a coefficient —2fl ^M& — e -^p 2 u ^ we ca jj ^ an internal 2n-leg cosine vertex) and 
the external or internal 2n-leg vertex. We will graphically represent the internal 2n-leg 
cosine vertex by a dot-centered circle with In legs. 

From the analysis and results in last paragraph, the sum of all e2 n -C2n graphs can be 
depicted in Fig. 3. In Fig.3, r = and a 2n-hne graph has a factor (2n)!. Note that in the 

In lines 

+ (x©® + ®Q® + ••• + +••• = ®o® 

rb ri ro r\ ro fi fo ri fo fi 

FIG. 3: The sum of all connected first-order graphs consisting of e2 n -vertices and the C2 n -vertices. 
right hand side of Fig.3, we have used the graph with two bold lines to represent the sum 
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of all possible graphs with even lines between two vertices. Thus, the sum in Fig. 3, G^\ T , 
can be calculated as 

' 2n [iV8^(3) 2n 



G% = Xle^G® / d»n £ 



7ra 



■(/« 



l\2n ' 
Ori . 



n=l 



(2n)! WOri ' (2n)! 
2ne- 4n{a2+(32)u J d u n[cosh(87ra/3f-^) - 1] . 



(2n) 



(47) 



Comparing last equation and the right hand side of Fig. 3, one see that the two bold lines 
corresponds to the integrand [cosh(87ra/3/ ^.|) — 1] and the factor before the integral in Eq.(47) 
comes from the n-independent parts of the coefficients for the external 2n-leg exponential 
and internal 2n-leg cosine vertices. It is evident that summing up in Eq. (46) and G^ n 
in Eq.(47) gives rise to Eq.(35). 

At the second order, diagrams will be complicated. There exist nine types of graphs: 
two types for one e2 n vertex connecting with two //-vertices (figures (a) and (b) in Fig.4), 
four types for one e 2n vertex connecting with two C 2n - vertices (figures (f), (g), (h) and (i) 
in Fig.4), and three types for connections among one C2 n vertex, one /i- vertex and one Ci n - 
vertex (figures (c), (d) and (e) in Fig.4). The sum of every type of graphs is depicted in 
Fig.4. In Fig.4, the graph with three bold lines between two vertices represents the sum of 
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FIG. 4: The sum of all connected second-order graphs consisting of e2 n -, C2n-vertices and //-vertex. 



all similar graphs with odd lines between the same two vertices, and a set of three bold lines 
corresponds to a negative hyperbolic sine function (one can check it by doing as was done 
in Eq.(47)). The results of graphs in Fig.4 can be calculated as follow, respectively. The 



27 



graph (a) in Fig. 4 is simply 

1 (i^aY 



a.a 



2 4! 



= 8vr 2 aVe- WW | dT? x j d"r 2 (/^) 2 (/^) 2 (48) 

In Eq.(48), the symbol C| is a number of combinations, and we have added the factor \ 
which is the total factor at the second order. We also do so for the other second-order 
graphs. The graph (b), (c), (d), (e), (/), (g), (h) and (i) in Fig.4 can be written as 

\2 



G% = Y^^Gf J d»n J ^ 2 /^^(~^ 2 )(~^)-^-2-2 

= - 47 raV 4 e- Ww / dTP, J d^f^fzl , (49) 

= \^^Gf J cfrx y rfV 2/ - /r - {-sinh[ 87 ra/3/-i]} 

x(-^ 2 )(-2Q) ( ^ )1 e-^-2-2 
= 16naf3f, 2 Qe~^ a2+ ^ u J d»n J d^r^jZ^si^aftf^] , (50) 

GS = d v r x J ^r 2 [cosh(87ra/3/-i) - 1] 

x(-2Q)i^^e- 4 ^(/- 2 ) 2 (-i^) -2-2 
= -87r/3y^e-M a2+ ^ w | ^ | d"f?,(/^) 2 [cosh (8™/%; ) - 1] , (51) 

Gg) = i | d"ft | rfV 2 (-2fi)e- 4 ^[cosh(87ra/?/ - r .;)-l] 



xG (0)(W (/ - )V l, 2) . 2 . 2 



= -St^V^- 4 ^^ | d"fi | rfV 2 (/-4) 2 [cosh(87ra/3/ >;) - 1] , (52) 

= J d»n J rfV 2 (-2fi)e- 4 ^[-sinh(8vra/5/-,;)] 

x(-2fi)e- 4 ^[-sinh(87ra^)]hsinh(8vr/3 2 / r ^ 2 )] 
= -2Q 2 e- 47r(a2+2 ^ u J d»n J d v r 2 smh(8na[3f^) 

x sinh(87ra/3/ >;) smli^ 2 /^) , (53) 
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Gg> = \g® j d»n J rfV 2 (-2fi) e - 4 ^[cosh(87ra/?/ >;) - 1] 

x(-2fi)e- 4 ^ 2w [cosh(87ra/3/- r ,;) - 1] [cosh(87r/5 2 ^i ) - 1] 
= 2Q 2 e- 4 ^ a2+2 ^ u J d u nj tff 2 [cosh(87ra/3/^) - 1] 

x[cosh(87ra/3/-i) - l][cosh(87r/? 2 /:: 4 2 ) - 1] , (54) 

G { ai = \G® Jd»n I rfV 2 (-2fi)e- 4 ^[cosh(8 7 ra/3/ ;i)-l] 

x(-2fi)e- 47r ' 32w [cosh(87r/? 2 /r 1 , 2 ) - 1] 
= 2Q 2 e- 4 ^ a2+2 ^ u J d u nj tff 2 [cosh(87ra/3/^) - 1] 

x[cosh(8vr/3 2 /,7 4 2 )-l], (55) 

GS = \j d»n J ^r 2 (-2fi)e- 4 ^ 2w [cosh(87ra/3/--.;) - 1] 



x 



Gi°)(-2fi)e- 4 ^ 2w [cosh(87ra/3/-i) - 1] 



= 2fi 2 e- 47r ( a2+2/32 ^{ J tffi[cosh(87ra/3/^) - l]} 2 . (56) 

For easily reading, every first expression in Eqs.(48) — (56) was written down in the way of 
one factor in the expression corresponding to one part in the relevant graph (some graphs 
have symmetrical factors). One can check that the sum of Eqs.(48) — (56) coincides with 
Eq.(39). Thus, one has seen that the VP expansion of one-point functions can be performed 
by borrowing Feynman diagrammatic technique in conventional perturbative theory. 



IV. SG AND SHG FIELD THEORIES WITH v = 2: ONE-ORDER VP RESULTS 

Subsection A and B will provide comparisons of the one-order VP results for sG and shG 
field theories with the conjectured exact results, respectively, and subsection C will compare 
the one-order VP results with the one-order perturbative results sG field theories . 



A. Comparisons with the Conjectured Results: sG Field Theory 

In two-dimensional Euclidean space, Eqs.(32) and (35) give (e = 1) 

n=l ^ ' 
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^ W 'l~2a 2 K 02 )+4n^(^f^K 0c , (57) 



y M 2J v ' M 2X M 2 

where, 



roo roc 

K 02 = / cIxxKq(x), K 0c = / <ixx[cosh(4a/3i^o( a; )) — !]■ 
Jo Jo 



To determine the auxiliary parameter /x, as stated in Sect. II, we can require that = 
according to the PMS, and have 

u 2 / Q l-a 2 -3 2 XVCI-/* 2 ) 

' f ^Koc -± \Z J • (58) 



7W 2 V a 2 (l -2a 2 K 02 ;/ 

Thus, Eq.(57) with Eq.(58) gives G* a , the approximate result of G a up to the first order in 
the VP theory. Note that the normal-ordering parameter Ai can be taken as any positive 
value with mass dimension, and different choice of its value leads to just a finite multi- 
plicative redefinition of Q 25] (1975). Once At is chosen, the theory is defined and the 
parameter Q is given a precise meaning, just as the choice of the normalization of the field 
cos(/?0) in Ref. |4| does. We will take Ai as unit. This choice renders every quantities in 
Eqs.(57) and (58) dimensionless, and simultaneously amounts to taking the same normaliza- 
tion conditions as Eqs.(6) and (16) in Ref. Q. The latter point can be checked by directly 
calculating the relevant two-point correlation function with the two points approaching each 
other within our formalism and then comparing the result with that in Ref. j4j. Thus, our 
results can be compared with the exact formula. Because being exponential-like interaction, 
the two-dimensional sG field theory can be rendered finite only for the case [3 2 < I by the 
Coleman's normal-ordering prescription (it amounts to the renormalization of the mass pa- 
rameter. For the range of (3 where the sG field theory diverges, one has to resort to a further 
renormalization procedure), and accordingly G a is finite for the case (3 2 < | and a/3 < \. 
And so G^ here is finite only for the range a/3 < | (one can see it by noticing the integral 
Kq c ). Fortunately, this range of (3 is just the validity scope of the conjectured exact formula 
in Ref. Q. Numerically, G ! a can be given with Mathematica programme for the range of 
a[3 < 0.426925, as stated in Ref. This also occurs at the shG field theory. 

Now we consider the comparison for the sG field theory. From Eq.(12) in Ref. jj], the 
dimensionless Q can be written as 



n 



Y{(3 2 ) 



tiY(1-P 2 ) 



2r(|) 



2-20- 



■2 



(59) 
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with £ = j-rg2- 111 l as t equation, we took the soliton mass as unit to make f2 dimensionless. 
For the conjectured exact formula Eq.(20) in Ref. 4], we also do so and have 

2sinr^)m + $)r(i-|) 



G a 



exact _ V 2 / \2 1 2' V 2' 



4^ 



r [°° dt sinh 2 (2a/3t) 2-2*1-1 

eXP 7o T [ 2sinh(^)sinh(t)cosh((l-/5 2 )t) " 2a 6 J} " (60) 

To compare with the conjectured exact result, we depict according to Eq.(57) with 

Eqs.(58) and (59) as the left figure of Fig. 5, and the exact result Gf act according to Eq.(60) 

as the right figure of Fig.5. In Fig.5, the parameter range is < a < 0.65 and < (3 < 0.65. 




FIG. 5: Comparison of G J a (the left) with the conjectured exact result G exact (the right). 

In principle, the two figures in Fig.5 can be extended to all range a/3 < 0.5. From Fig.5, one 
can see that for smaller f3 and larger a, G a tends to zero, and for larger (3 and a, G a is greater 
than the value 1. When a large f3 is given, G a increases with the increase of a, and when a 
small j3 is given, G a decreases with the increase of a. On the other hand, the two figures in 
Fig.5 resemble each other very well and suggest that for the range of the parameter {a, /3} 
drawn in Fig.5 G T a has a good agreement with the conjectured exact result. Our numerical 
analysis indicates that for the case of a < 0.2 or so, the one-order results almost completely 
agree with the conjectured exact results, and for larger a, our results differ from the exact 
results with about ten percents or so (at most with 20 more percents when a approaches 
values with a[3 < 0.426925 satisfied, see Table I in Ref. For a concrete illustration, 

we draw G ! a (the solid curves) and the exact result G^ xact (the dotted curves) at (3 — 0.15 
and 0.7 in Fig. 6 and Fig. 7, respectively (the dashed curves is the perturbative results, see 
next subsection C). These two figures indicate that G T a is always greater than G exact . 
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FIG. 6: Comparison between the one-order (G^, solid curve) and conjectured exact (G^ xact , dotted 
curve) sG VEVs of exponential fields at (3 = 0.15. The dashed curve is the one-order perturbative 
result (G I a pert , see subsection C) at the same value of (3. At this value of j3,G I d pert and G^ xact are 
widely discrepant. 




0.2 0.4 0.6 

FIG. 7: Similar to Fig.6 but at (3 = 0.7. At this case, the one-order perturbative result (the 
dashed curve) is almost coincident with the conjectured exact result. 

B. Comparisons with the Conjectured Results: shG Field Theory 

In the Euclidean space, the shG Lagrangian density is 

£shG — -V 'r4>r + 2A COsh(v / 87T70r) • (61) 

2 
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For the convenience of comparison with the conjectured exact VEV of the exponential field 
in the shG field theory in Ref. jfj], here we are interested in the VEV Gb =< e^* b ^ >. 
Taking a = —ib, (3 = 27 and Q = —A in Eq.(57), we get the one-order VP approximate 



result of Gb, G{, as follows 



G{ = ^ 2 Y b \l + 2b 2 K 02 ) - 4nA( f i 2 )- b2 ^ 2 - 1 K 0ch 



(62) 



with the parameter fi 2 satisfied 



U 



(63) 



where, K 0c h = / °° dxx[cosh(4:b'yKo(x)) — 1] is similar to K Qc . From Eqs.(9) and (8) in 
Ref. |5[, the dimensionless A has the following form 



A 



r(- 7 2 



7rr(l + 7 2 ) _ 



'(2+2 7 2 ) 



)rn + 



2+2 7 2 



4V? 



2+2 7 2 



(64) 



and the dimensionless conjectured exact expression of Gb, G e h xact , is 



r^exact 



1 \(2+2-y 2 )) L \ L ^ 2+2 7 2 



exp 



00 dtr 



sinh 2 (2&7t) 



2b z e 



2 -2t 



]}. (65) 



,„ t L 2 sinh(7 2 t) sinh(t) cosh((l + 7 2 )t) 
In Eqs.(64) and (65), we adopted symbols here and took the particle mass as unit. Now we 
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FIG. 8: Comparison of G T b with the conjectured exact Gb- 

can compare G^ with G e h xact . We depict G{ (the left figure) and G e h xact (the right figure) in 
Fig. 8 for the range of b < 0.9 and 7 < 0.9 (it can be extended to the all tractable range of 
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67 < 0.426925 ). In the left figure of Fig.8, G{ is set to zero for the range of 67 > 0.426925 
where Mathematica couldn't produce finite numerical results. The two surfaces in Fig.8 are 
basically resemble each other, except for the cases of larger b and simultaneously smaller 
7. Numerical analysis indicates that for a given value of j3, G[ always decreases with the 
increase of b, while G e h xact is first increase and then goes down with the increase of b. But, 
for not large b, the relative differences of G[ from the conjectured exact result Gl xact are 
small, the situations on the differences between G{ and Gl xact are analogous to those for the 
sG field theory. We show these points in Fig. 9 and Fig. 10 for the values of 7 = 0.2 and 




0.5 



1.5 



FIG. 9: Comparison of G\ (solid) with the conjectured exact G& (dotted) at 7 = 0.2. 
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FIG. 10: Comparison of G\ (solid) with the conjectured exact Gb (dotted) at 7 = 0.8. 



7 = 0.8, respectively. In Figs. 9 and 10, the dotted curves are for G% and the solid ones 
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for G b . Specially, when b is smaller than 0.2 or so, G\ almost completely coincides with the 
conjectured result G e h xact . These two figures show that for the shG field theory, G{ is always 
smaller than G e h xact , which is opposite to that for the sG field theory. 

By the way, if we choose A as unit instead of setting the particle mass as unit as did 
in the above, then the numerical analysis indicates that the dependence of G[ upon b at a 
given (3 is qualitatively similar to that of G e h xact . 

C. Comparisons with the One-order Perturbative Results: sG Field Theory 

Using the massive Thirring model, Ref. [3] (Poghossian) calculated G a perturbatively up 
to the first order of the coupling g, and obtained 



= (> 2 exp{ r^[ Sinh2( f at) -2aV^]} 
V P Vo t 1 sinh 2 (t) JJ 



x [1 + A [_ 2 tf (I) + 8a 2 (l - 4 log 2) + 8a 2 (^( v / 2a) + *(->/2a)) 

07T 2 

+ (1 - 8a 2 )(*( ^— ) + *( ^— ))]] + o(<? 2 ) , (66) 

where, the coupling g is related to the sG coupling (3 as follows 

l = (67) 

Eq.(67) indicates that g — > corresponds to /3 2 — f ^. In Fig. 11, we compare the one-order 




0.2 0.4 a 6 



FIG. 11: Similar to Fig.6 and Fig.7, but at (3 = 0.65. 
perturbative (the dashed curves) and VP (solid) results with the conjectured exact results 
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(dotted) at the case of (3 — 0.65. In subsection A of this section, the same comparisons 
were made at (3 — 0.15 and 0.7 in Fig. 6 and Fig. 7, respectively. In Figs. 7, (3 2 = 0.49, 
i.e., (3 approaches \ or g — > 0, the one-order perturbative results are almost completely 
identical to the conjectured exact results for all the plotted values of a, while the one- 
order VP results have evident differences (the relative errors are 20% or so at most) from 
the conjectured exact results for larger values of a. In Fig. 11, f3 2 ~ 0.42, the one-order 
perturbative results (dashed) have quite large deviations from the conjectured exact results 
(dotted) for larger values of a, and differences of the one-order VP results (solid) from the 
conjectured exact results for larger values of a get smaller than the case of (3 2 = 0.49. In 
Fig. 6, f3 2 = 0.02, the one-order perturbative results and the conjectured exact results are 
widely discrepant for all the plotted values of a, while differences of the one-order VP results 
from the conjectured exact results for larger values of a get much smaller than the cases of 
f3 2 = 0.49 and 0.42. Besides, in Figs. 6, 7, and 11, the one-order VP results almost coincide 
with the conjectured exact results for not large a (< 0.2 or so). To illustrate the dependency 
of the deviations of the one-order perturbative and VP results from the conjectured exact 
results upon the coupling f3, G Ipert , and G e * act with a = 0.2, 0.3 and 0.5 are depicted 
as functions of (3 in Figs. 12, 13 and 14, respectively. These three figures evidently show 




FIG. 12: The dependances of G{ pert (dashed), G^(solid) and G e a xact (dotted) upon f3 for a = 0.2. 

that the one-order perturbative results (dashed) have good agreements with the conjectured 
exact results (dotted) only when j3 2 approaches the value | (i.e., when g — > 0), but their 
deviations from the conjectured exact results appear, for a given value of a (not too small), 
and become larger and larger when (3 2 decreases from \. For a given a, when values of (3 2 
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G a 




0.2 0.4 0.6 

FIG. 14: The dependances of G 1 ^ 1 (dashed), G{ (solid) and G e a xact (dotted) upon (3 for a = 0.5. 



is much smaller than |, the one-order perturbative and conjectured exact results are widely 
discrepant. On the other hand, the one-order VP results (solid) are almost identical to the 
conjectured exact results for not large values of (3, and have not too larger differences from 
the conjectured exact results even when f3 2 approaches \. 

V. CONCLUSION 

In this paper, we developed a VP scheme to calculate the VEVs of local fields in relativis- 
tical QFTs. For a class of scalar field theories whose potential have Fourier representations, 
we obtained the Gaussian smearing formulae for the VEVs of a comparatively general local 
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field. As an application and illustrations on the scheme and the Gaussian smearing formu- 
lae, we considered the sG and shG field theories. This application provided an example in 
both directly and Feynman-diagrammatically performing the VP scheme, and showed the 
usefulness of the Gaussian smearing formulae. The Gaussian smearing formulae can relieve 
us of hard labor in path integrals. The usefulness of the Gaussian smearing formulae can 
also be envisioned from the smearing formulae for the higher order effective classical po- 
tential in statistical mechanics {^(J, which succeeded in applying to the singular Coulomb 
potential. Additionally, according to the saying in Ref. j^J when discussing normal-ordering 
prescription, the existence of the Fourier representation of a given V(4>) for the validity of 
the Gaussian smearing formulae is irrelevant in the derivation of section II, because the final 
Gaussian smearing formulae are purely algebraic [2]. Although the one-order VP results 
are the lowest approximate in VP theory, the numerical discussions in last section indicated 
that the one-order VP VEVs of the exponential fields for the sG and shG field theories give 
a strong support to the conjectured formulae on the exact VEVs in Refs. and the 

first paper in Ref. 0, at least, for not large values of the coupling and the exponential-field 
parameter. They also suggest the effectiveness of the VP scheme here in calculating the 
VEVs of local fields for QFTs. The comparisons in last section also illustrates the non- 
perturbability of the VPT scheme here and its advantages over the perturbative theory. We 
believe that the VP scheme in the present paper can provide an effective, systematical con- 
trollable non-perturbative approximate tool for calculating VEVs of local fields in relativistic 
QFTs. 

As was pointed out in Ref. j^, there exist some interesting problems based on our 
work, and here we do not repeat them one by one. Nevertheless, we intend to stress that 
a further investigation on the higher order results would improve the one-order VP results 
here and show the convergency of the VP theory. As a matter of fact, a numerical analysis 
on the VP results up to the second order in the sG field theory in Ref. [3] has suggested 
this point. The applications of the scheme in the present paper to other field theories 
maybe also give a good check and substantial support to the other existed conjectured exact 
formulae Simultaneously, the problem of calculating VEVs of 

local fields provides a good laboratory for the VP theory because there existed so many 
conjectures on VEVs of local fields in various perturbed conformal field theories. Of course, 
generalizing the VP scheme here to other physical problems will be interesting and useful. 
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For example, the bound state problems in relativistic QFT is a notoriously difficult non- 
perturbative problems, and we noticed that there had existed an approach of attacking it 
which is a combination of variational method and perturbative theory j^]. The approach in 
Ref. can be regarded as a method of calculating VEV, but, as pointed by the authors, 
it is valid only for the weak coupling, because the variational procedure in it is to choose 
the appropriate two-particle operator and the expansion scheme is a naive perturbative 
expansion. We think that based on the work in this reference, the VP scheme here can be 
generalized for calculating bound-state mass, and this generalization will be meaningful and 
interesting. 
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